Quantum simulation of a system with competing two- and three-body interaction^] 
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Quantum phase transitions occur at zero temperature, when the ground state of a Hamiltonian 
undergoes a qualitative change as a function of a control parameter. We consider a particularly 
interesting system with competing one-, two- and three-body interactions. Depending on the relative 
strength of these interactions, the ground state of the system can be a product state, or it can 
exhibit genuine tripartite entanglement. We experimentally simulate such a system in an NMR 
quantum simulator and observe the different ground states. By adiabatically changing the strength 
of one coupling constant, we push the system from one ground state to a qualitatively different 
ground state. We show that these ground states can be distinguished and the transitions between 
them observed by measuring correlations between the spins or the expectation values of suitable 
entanglement witnesses. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Mn 



Introduction.- At zero temperature, a system can un- 
dergo a quantum phase transition (QPT) to a new 
ground state as a result of a change in a param- 
eter of the Hamiltonian [1 j. Well-known examples 
are the superconductor-insulator transition and the 
paramagnetic-antiferromagnetic transition in quantum 
magnets. QPTs were experimentally observed in mag- 
netic systems [2] , heavy- fermion metals [3 j , common met- 
als [4] and Bose-Einstein condensation [5 j. The investi- 
gation of QPTs is useful for discovering novel materials 
[6], and some QPTs show interesting entanglement char- 
acteristics [7 j. A collection of reviews [6, 8 on the topic 
of QPTs reported the current status and recent develop- 
ments in this field. 

In most systems studied, attention was focused on 
two-body interactions, which are most readily accessi- 
ble experimentally. On the other hand, systems with 
three-body interactions have been shown to exhibit ex- 
otic quantum phases in their ground states [9j, such as 
topological phases or spin liquids and a chiral phase. 
However, it is difficult to observe these properties in ex- 
periments. The main challenges are: (i) to identify ex- 
perimentally accessible systems with three-body interac- 
tions; (ii) to experimentally control the variation of the 
system Hamiltonian in a sufficiently precise manner; (Hi) 
to characterize the resulting ground state. 

QPTs are generally associated with strongly correlated 
quantum systems, which cannot be efficiently simulated 
on classical computers because the required computa- 
tional resources grow exponentially with the system size. 
As suggested by Feynman [10] and proved by Lloyd [IT], 
however, a quantum computer can efficiently perform this 
kind of simulations and provide new insight into strongly 
correlated quantum systems, including QPTs. Previ- 
ous examples of such studies include the simulation of 
a three-body Hamiltonian [12 and quantum magnets in 
a two-spin Heisenberg chain using nuclear spins [13] as 



well as trapped ions [T4] . 

In this Letter, we experimentally simulate the small- 
est spin system involving a three-body interaction in an 
NMR quantum simulator. Such systems are the building 
blocks of triangular Ising nets [15] and triangular spin 
ladders [16] . which have been studied in detail in the 
context of quantum statistics [17] and condensed matter 
theory [9] [15]. In general, the competition between the 
different interactions results in QPTs. Here we prepare 
the system in the ground state of the Hamiltonian and 
drive it from one phase into a different phase by adiabati- 
cally changing the Hamiltonian. By quantifying different 
types of entanglement, or by using suitable entanglement 
witnesses [19|, we successfully detect the quantum tran- 
sition induced by the three-body interactions. In the 
thermodynamic limit, this transition should correspond 
to a novel type of QPT [T6] , 

System. - Consider a chain of N spins 1 /2 in a uniform 
magnetic field, interacting by Ising-type nearest-neighbor 
two-body and three-body couplings: 

+J 2 ^^ 1 + J 3 ^^V^, (1) 

where the o l z are the Pauli operators, uo z and u x the 
strengths of the longitudinal and transverse magnetic 
fields, and J2, J3 the two-body and three-body coupling 
constants. The competition between the three kinds of 
interactions (one-, two- and three-body) determines the 
ground state of the system. We first consider two limit- 
ing cases, corresponding to the two- and three-spin Ising 
models. In both cases, we discuss explicitly the situation 
for positive coupling constants (jj z ,uj x , J2, J3 > 0. 

(z) Two-spin Ising model (J3 = 0) [1 j: {id) in a lon- 
gitudinal magnetic field (^=0). The level crossing at 
J2 = uo z corresponds to a first-order phase transition of 
the ground state from a paramagnetic state \ijjg) = !•• ii 
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..) in the weak-coupling case (J2 < uj z ) to a two- fold 
degenerate, antiferromagnetically ordered ground state: 
\^f F ) = {|.. U ••),!•• It ••)} in the strong-coupling 
case (J2 > uo z ). In the case of ferromagnetic coupling 
(J2 < 0), no QPT occurs, (ib) Two-spin Ising model in a 
transverse magnetic field (i.e., uj z =0). This well-studied 
model exhibits a second-order QPT at J2 = uo x in the 
thermodynamic limit pQ , where the ground state changes 
from the paramagnetic phase \ip P ) = |.. ^— ^— ..) to the 
antiferromagnetically ordered, doubly degenerate ground 
state \^j F ). Here | <-) = (| t) - | l))/V2. These two 
cases are exactly solvable and the induced phase transi- 
tions can be detected by the traditional two-point corre- 
lation functions pQ. 

(ii) Three-spin Ising model with J 2 = in a transverse 
magnetic fields (i.e., cj z =0): This model is not exactly 
solvable, but numerical simulations predict a critical 
point at J3 = uj x . The strong-coupling phase ( J3 > uj x ) is 
fourfold degenerate: \ijj F ) = {|.. ||| .), |.. J/ft ..), |.. tit 
..), |.. ttt ••)}? while the non-degnerate ground state of 
the weakly-coupled phase is the paramagnetic state \ijj P ) 

In the following, we demonstrate these features in a 
proof-of-principle experiment using N = 3 spins and peri- 
odic boundary conditions cr^ +1 = o\. Figure [l] shows the 
ground state of the system. The competition between the 
three kinds of interactions (one-, two- and three-body) re- 
sults in different ground state phases (product-, W-, and 
GHZ-states). In regime I, the ground states are prod- 
uct states, in phase II, they are W-type states, and in 
regime III, they are GHZ-type states. The GHZ and W 
states are the only two inequivalent kinds of genuine tri- 
partite entanglement in a 3-spin system [21]. As a result, 
such a simple model allows us to experimentally study 
the two inequivalent states of genuine tripartite entan- 
glement and the novel phenomena induced by three-spin 
interactions. 

Quantum simulation.- We chose the Diethyl- 
fluoromalonate molecule as a three-spin NMR quantum 
simulator, in which the 13 C, X H, and 19 F nuclear spins 
are represented by spins 1, 2, and 3, respectively. The 
molecular structure and the relevant parameters are 
shown in Fig. [2] (a). The natural Hamiltonian of the 

system is U NM r = ELi + Ei <jM =i ^^M- 

From this natural Hamiltonian, we generate the model 
Hamiltonian ([!]) with equal two-body coupling strengths 
by a suitable refocusing scheme (22j. The three-body 
interaction can be simulated by a combination of two- 
body interactions and RF pulses [12] . Since all terms 
in the Hamiltonian ([!]), except the transverse field 
term, commute with each other, we expand the over- 
all evolution by the following concatenation: e~ inr = 

e -m x r/2 e -m z r e -m x r/2 + Q^y ^ = £ <T» 

and 1~L Z = T~L — H x . This expansion faithfully represents 
the targeted evolution provided the duration r is kept 
sufficiently short. Fig. [2^b) shows the pulse sequence 
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FIG. 1: (Color Online) Schematic phase diagram of the 
ground state of the Hamiltonian % for N — 3. Here 
|GHZ±) = (| ttt) ± I iU))/y/2 and H is the Hardamard 
gate. The arrows in (a) and (b) represent the adiabatic evo- 
lutions discussed in the experimental section. For (a) and (c), 
uj x is assumed to be small, i.e., <^ \cj z \. 
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FIG. 2: (Color Online) (a) Molecular structure and proper- 
ties of the quantum register: Diethyl-fluoromalonate. The 
oval marks the three spins used in the experiment. The ta- 
ble on the right summarizes the relevant NMR parameters 
measured at room temperature on a Bruker Avance II 500 
MHz (11.7 Tesla) spectrometer, i.e., the Larmor frequencies 
uji/2n (on the diagonal), the J-coupling constants Jij (below 
the diagonal), and the relaxation times Ti and T2 in the last 
two columns, (b) Pulse sequence for simulating the Hamil- 
tonian of Eq. (TlJ) . The narrow black rectangles represent 
small-angle rotations, the narrow empty rectangles denote 
90° rotations and the wide ones denote the refocusing 180° 
pulses. The delays are n = J2T /[l/ (nJij) + l/(jrJjk)] with 
(i, j, k) an even permutation of (1, 2, 3) and d\ — 2 J37-/ (7r J12). 
The offsets between the irradiation frequencies and the corre- 
sponding Larmor frequencies are FQ1 = 2uj z t/(t\ — 72 + 373), 
FQ2 = 2uj z t/{t 1 + 7-2 - 73) and FQ3 = 2u z t/(ti + t 2 + r 3 ). 



that realizes e~ lUr . 

In order to observe the system undergoing the tran- 
sitions, it must always be close to the instantaneous 
ground state of the time-dependent Hamiltonian. This 
was achieved by quantum adiabatic evolution [23] , which 
requires: (i) the system starts in the ground state of 
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H(0), (ii) H(t) changes sufficiently slowly to satisfy the 
adiabatic condition. For the experimental implemen- 
tation, we discretized the time-dependent Hamiltonian 
H(t) into M + 1 segments H(m) = n[C(§T)} with 
m = 0,1,..., M [24], where T is the total duration of 
the adiabatic passage and C is the control parameter 
in H(t). The adiabatic condition is satisfied when both 
T, M — ^ oo and the duration of each step r — >• 0. 

To observe the different types of transitions in the 
Hamiltonian 0, we chose two different parameter sets: 

Case A: Transition from a product state to a W-type 
entangled state: the control parameter is the 
two-body coupling strength J2, which varies from 
^2(0) = to J2CO = 2. The other parameters of 
the Hamiltonian are constant, uj z = —2, uj x = 0.09 
and J3 = 0, as shown in Fig. [I] (a). 

Case B: Transition from a product state to a GHZ- 
type entangled state: the control parameter is the 
three-body coupling strength J3, which varies from 
J 3 (0) = to J 3 (T) 2, with uo z = J 2 = and 
uj x = 0.12, as shown in Fig. [I] (b). 

We used a hyperbolic sine for the time dependence of the 
control parameters J2W and J%(t) in the experiments 

To determine the optimal number M of steps in the 
adiabatic transfer, we used a numerical simulation of the 
minimum fidelity encountered during the scan as a func- 
tion of the number of steps into which the evolution is 
divided (see Fig. [3|. The fidelity is calculated as the 
overlap of the state with the ground state at the relevant 
position. Using the decoherence model of Vandersypen 
et al [25], we also simulated the effects of decoherence 
under the actual (piecewise constant) Hamiltonian gen- 
erated by the pulse sequence of Fig. [5](b). The resulting 
fidelity of these points, which is represented by the boxes 
in Fig. [3J is consistently lower than that of the ideal 
Hamiltonian and reaches a maximum for M w 8. For 
the experiments, we therefore chose M = 8. 

Detection of transition points.- The conventional ap- 
proach to detect the system undergoing an entangling 
phase transition is to measure two-spin correlations like 
C xx — \ s ^i^j(°~x°~i) • The result of such a measure- 
ment is shown in Fig. [4] The result clearly shows the 
expected transition for Case A: The two-spin correlation 
goes through a step-like increase at J 2 ~ 1, consistent 
with the expectation that the ground state should change 
from a product state to a W state at this value of the cou- 
plings strength. However, in Case B, where the system 
should start in a product state and end up in a GHZ 
state, the measured two-spin correlations give no indica- 
tion of this transition. 

In order to obtain a clean signature of this novel tran- 
sition, induced by three-body interactions, we employed 
entanglement witnesses [19] . These witness operators can 
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FIG. 3: (Color Online) Numerical simulation of the minimum 
fidelities during the adiabatic passage vs. the number of steps 
for the ideal (piecewise constant) Hamiltonian (O) and for the 
more realistic case that includes the actual pulse sequence and 
the effect of decoherence (□) for both cases A and B. 
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FIG. 4: (Color Online) Experimental detection of quantum 
transition points by two-spin correlations C|^ p (* and +), to- 
gether with the simulated correlations C^™ (O an d □), which 
take the effect of decoherence into account. The effective de- 
coherence time was estimated as 150 ms for Case A and 
600 ms for Case B. 



always be used to detect various forms of multipartite en- 
tanglement, provided we have some a priori knowledge 
about the states under investigation. In our context, the 



and 



two witness operators Ww = §1 — |Wooi)(Wooi| 
Wqhz = fl - |GHZ_)(GHZ_| are useful for detecting 
multipartite entanglement. If their expectation value is 
negative, they indicate the presence of genuine tripartite 
entanglement; Tr(p>VGHz) < specifically detects GHZ- 
type entanglement. To perform these measurements, we 
applied a basis transformation to the density operator 
and destroyed the off-diagonal elements with pulsed mag- 
netic field gradients to implement the projective measure- 
ment [27 . The populations were read out by applying, in 
three separate experiments, tt/2 readout pulses to each 
qubit and measuring the resulting transverse magnetiza- 
tion. 

We first initialized the system into the ground state of 
the initial Hamiltonians [26] and then varied the Hamil- 
tonians adiabatically. To minimize nonadiabatic effects, 
the scan had to be relatively slow, lasting about 146 ms 
for Case A and 62 ms for Case B. Spin relaxation there- 
fore reduced the signal for the later times compared to 
the initial part of the scan. To separate this effect from 



4 



the effect of the transition and obtain a cleaner signature, 
we first determined the overall signal decay and rescaled 
the data to remove this effect. 

Fig. [5] summarizes the experimental results for both 
cases, using both entanglement witnesses. Both transi- 
tions are now clearly visible, and the tripartite entan- 
glement is present in both final states. Together with 
the measurement of the two-spin correlations in Fig. |4j 
we can also differentiate between the two types of transi- 
tions. In Case A, the increasing strength of the two-body 
interaction J2 pushes the system from a product state to 
a W state, but in Case B, the increase in the three-body 
interaction strength J3 generates a GHZ state. 
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FIG. 5: (Color Online) Experimental detection of quantum 
transition points by entanglement witnesses for Case A (left 
hand side) and Case B (right hand side). The solid lines 
represent the theoretical expectations. 

We also determined the details of the final state at 
the end of the scan, by performing complete quantum 
state tomography [28]. The experimentally generated 
state p exp had an overlap of F = \(^id\Pex P \^id)\ with 
the ideal state \ijjid) of 0.61 for the W state and 0.73 
for the GHZ state. The reduced fidelity mainly results 
from relaxation during the long adiabatic passage. As in 
the case of the witness measurements during the scan, 
we also removed the effect of spin relaxation from the 
measured final states by calculating the 'experimental fi- 
delity' [29] F = \{^ id \p exv \^id)/Tr(pl xp )\. The observed 
data yielded F(p™ p ) = 0.90 and F(p™ z ) = 0.92, in 
excellent agreement with theoretical expectations. 

Conclusion.- Using an NMR quantum simulator, we 
have performed an experimental quantum simulation of 
a system with competing one-, two- and three-body 
interactions. By adiabatically changing the Hamilto- 
nian, we could observe the system undergoing two dif- 
ferent entangling transitions. The generation of entan- 
glement confirms that the transitions are of quantum- 
mechanical origin (analogous to quantum fluctuations in 
QPTs). In particular, with suitable entanglement wit- 
nesses as observables, we detected a novel transition in- 
duced by three-body interactions that is qualitatively dif- 
ferent from the states that can be characterized by two- 
spin correlations. 



We consider our present experiment based on a uni- 
versal simulator as a first experimental step that demon- 
strates how one can obtain a clear picture of phase transi- 
tions in ground states of many-body Hamiltonians, with 
the help of tools coming from the theory of entanglement. 
Systems consisting of more than three subsystems can 
display even more different phases in their ground states. 
We expect that entanglement witnesses will turn out to 
be useful for detecting ground-state phases in larger sys- 
tems, together with other entanglement measures (e.g., 
the global entanglement [30]). 

We thank Y. J. Deng, D. Bruss and H. Kampermann 
for helpful discussions. This work is supported by the 
CAS, NNSFC, and the DFG through Su 192/19-1. 



* xhpeng@ustc.edu.cn, djf@ustc.edu.cn, Dieter. Suter@tu- 
dortmund.de. J. Zhang is presently at Institute for Quan- 
tum Computing, University of Waterloo, Canada. 
[1] S. Sachdev, Quantum Phase Transition (Cambridge Uni- 
versity Press, Cambrige, 1999). 
[2] H. M. R0nnow et al, Science 308, 389 (2005). 
[3] J. Custers et a/., Nature 424, 524 (2003). 
[4] A. Yeh et al, Nature 419, 459 (2002). 
[5] M. Greiner et a/., Nature 415, 39 (2002). 
[6] R C. Canfleld, Nature Phys. 4, 167 (2008); P. Gegenwart, 

Q. Si, and F. Steglich, ibid. 4, 186 (2008). 
[7] A. Osterloh et a/., Nature 416, 608 (2002). 
[8] Editorial, Nature Phys. 4, 157 (2008); T. Giamarchi 
et a/., ibid. 4, 198 (2008); S. Sachdev, ibid. 4, 173 (2008). 
[9] H. P. Buchler et a/., Nature Phys. 3, 726 (2007); 
C. D'Cruz and J. K. Pachos, Phys. Rev. A 72, 043608 
(2005); K. A. Penson, J. M. Debierre, and L. Turban, 
Phys. Rev. B 37, 7884 (1988); J. C. Angles d'Auriac and 
F. Igloi, Phys. Rev. E 58, 241 (1998); P. Suranyi, Phys. 
Rev. Lett. 37, 725 (1976); J. K. Pachos and E. Rico, 
Phys. Rev. A 70, 053620 (2004). 
[10] R. P. Feynman, Int. J. Theor. Phys. 21, 467 (1982). 
[11] S. Lloyd, Science 273, 1073 (1996). 
[12] C. H. Tseng et al, Phys. Rev. A 61, 012302 (1999). 
[13] X. Peng et al, Phys. Rev. A 71, 012307 (2005); J. Zhang 

et a/., Phys. Rev. Lett. 100, 100501 (2008). 
[14] A. Friedenauer et a/., Nature Phys. 4, 757 (2008). 
[15] G. H. Wannier, Phys. Rev. 79, 357 (1950). 
[16] J. K. Pachos and M. B. Plenio, Phys. Rev. Lett. 93, 
056402 (2004). 

[17] P. D. Francesco et a/., Conformal Field Theory 

(Springer, New York, 1997). 
[18] N. P. Ong and R. J. Cava, Science 305, 52 (2004); S. 

Nakatsuji et a/., Science 309, 1697 (2005); G. Aeppli and 

P. Chandra, Science 275, 177 (1997); M. Harris, Nature 

456, 886 (2008). 
[19] A. Sanpera et a/., Phys. Rev. A 63, 050301 (2001); 

M. Bourennane et al, Phys. Rev. Lett. 92, 087902 

(2004). 

[20] K. A. Penson et al, Phys. Rev. B 26, 6334 (1982). 
[21] W. Diir et a/., Phys. Rev. A 62, 062314 (2000). 
[22] N. Linden et al, Chem. Phys. Lett. 305, 28 (1999). 
[23] A. Messiah, Quantum Mechanics (Wiley, New York, 



5 



1976). 

[24] M. Steffen et a/., Phys. Rev. Lett. 90, 067903 (2003). 
[25] L. M. K. Vandersypen et al, Nature 414, 883 (2001). 
[26] D. G. Cory et a/., Physica D 120, 82 (1998). 
[27] D. G. Cory et a/., Phys. Rev. Lett. 81, 2152 (1998); 

G. Teklemariam et a/., ibid. 86, 5845 (2001);M. A. 

Nielsen et a/., Nature 396, 52 (1998). 



[28] I. L. Chuang et a/., Proc. R. Soc. Lond. A 454, 447 
(1998). 

[29] E. M. Fortunato et al, J. Chem. Phys. 116, 7599 (2002). 
[30] T. R. de Oliveira et a/., Phys. Rev. Lett. 97, 170401 
(2006). 



